BABAR-PUB-06/015 

SLAC-PUB-11799 

hep-ex/0604007 



Measurements of branching fractions, rate asymmetries, and angular distributions in 

the rare decays B — ► K£+t~ and B — ► K*l+l~ 
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We present measurements of the flavor-changing neutral current decays B — » K£ IT and B — > 
K*£ + £~ , where £ + £~ is either an e + e~ or /i + (i~ pair. The data sample comprises 229 x 10 6 Y(AS) — > 
BB decays collected with the BABAR detector at the PEP-II e + e~ storage ring. Flavor-changing 
neutral current decays are highly suppressed in the Standard Model and their predicted properties 
could be significantly modified by new physics at the electroweak scale. We measure the branching 
fractions B{B -»• K£ + £~) = (0.34±0.07±0.02) x 10~ 6 , B(B -> K'l+T) = (0.78+^;^ ±0.11) x 10~ 6 , 
the direct CP asymmetries of these decays, and the relative abundances of decays to electrons and 
muons. For two regions in £ + £~ mass, above and below mj/^ , we measure partial branching fractions 
and the forward-backward angular asymmetry of the lepton pair. In these same regions we also 
measure the K* longitudinal polarization in B — > K*£ + £~ decays. Upper limits are obtained for 
the lepton flavor- violating decays B — > Ke/j, and B — * K*efi. All measurements are consistent with 
Standard Model expectations. 

PACS numbers: 13.25.Hw, 13.20.He 



I. INTRODUCTION 

The decays B — * K^t + i~ , where £ + £~ is either an 
e + e~ or /i + /i~ pair and K^*> denotes either a kaon or 
the if* (892) meson, arc manifestations of b — > s£ + £~ 
flavor-changing neutral currents (FCNC). In the Stan- 
dard Model (SM), these decays are forbidden at tree level 
and can only occur at greatly suppressed rates through 
higher-order processes. At lowest order, three amplitudes 
contribute: (1) a photon penguin, (2) a Z penguin, and 
(3) a W + W~ box diagram (Figure ^). In all three, a vir- 
tual t quark contribution dominates, with secondary con- 
tributions from virtual c and u quarks. Within the Oper- 
ator Product Expansion (OPE) framework, these short- 
distance contributions are typically described in terms of 




FIG. 1: Examples of Standard Model Feynman diagrams for 
the decays B — > K^£ + £~. For the photon or Z penguin 
diagrams on the left, boson emission can occur on any of the 
b, t, c, u, s, or W lines. 
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the effective Wilson coefficients Cf, Cf, and Cjf Q. 
Since these decays proceed via weakly-interacting par- 
ticles with virtual energies near the electroweak scale, 
they provide a promising means to search for effects from 
new interactions entering with amplitudes comparable to 
those of the SM. Such effects are predicted in a wide va- 
riety of models HHSHQ. 

In the SM the B — > K£ + £~ branching fraction is pre- 
dicted to be roughly 0.4 x 10~ 6 , while the B -> K*£+£~ 
branching fraction is pre dicted to be about three times 
larger 0, H S E3, M d The B -> K*£+£~ mode 
receives a significant contribution from a pole in the pho- 
ton penguin amplitude at low values of q 2 = rnj +e _ , 
which is not present in B — > K£ + £~ decays. Due to the 
lower mass threshold for producing an e + e~ pair, this 
enhances the K *e + e" final state relative to the K* /z + pT 
state. Currently, theoretical predictions of the branching 
fractions have associated uncertainties of about 30% due 
to form factors that model the hadronic effects in the 
B — > K or B — <> K* transition. Previous experimental 
measurements of the branching fractions are consistent 
with the range of theoretical predictions, with experi- 
mental uncertainties comparable in size to the theoretical 
uncertainties jT^, [3 ■ 

With larger datasets, it becomes possible to measure 
ratios and asymmetries in the rates. These can typically 
be predicted more reliably than the total branching frac- 
tions. For example, the direct CP asymmetry 
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is expected to be vanishingly small in the SM, of order 
10" 4 in the B -> K*£ + £~ mode [H]. However it could 
be enhanced by new non-SM weak phases [lg. The ratio 
Rk, defined as 



R K 



T(B -> Kn+n~) 

T(B -> Ke+e-) ' 



also has a precise SM prediction of Rk = 1.0000 ± 
0.0001 [17| . In supersymmetric theories with a large ra- 
tio (tan/3) of vacuum expectation values of Higgs dou- 
blets, Rk can be significantly enhanced. This occurs via 
penguin diagrams in which the 7 or Z° is replaced with 
a neutral Higgs boson that preferentially couples to the 
heavier muons [lg. In B —> K*£ + £~ this ratio is mod- 
ified by the photon pole contribution, thus the SM pre- 
diction is Rk* ~ 0.75 |j| with an estimated uncertainty 
of 0.01 [I3 if the pole region is included, or Rk* ~ 1-0 if 
it is excluded |17j . 

Additional sensitivity to non-SM physics arises from 
the fact that B — ► K^£ + £~ transitions are three-body 
decays proceeding through three different electroweak 
penguin amplitudes, whose relative contributions vary 
as a function of q 2 . Measurements of partial branching 
fractions and angular distributions as a function of the 
invariant momentum transfer q 2 are therefore of particu- 
lar interest. The SM predicts a distinctive pattern in the 
forward-backward asymmetry 



Afb(s) 



L , n d 2 r (B^K { *h+r) . . 

dcostf - — Sign(cos( 



-1 



d cos 9 ds 



dT(B^ ifW £+£-)/ ds 



where s = q /nig, and 9 is the angle of the lepton with 
respect to the flight direction of the B meson, measured 
in the dilepton rest frame [lj|. In the presence of non-SM 
physics, the sign and magnitude of this asymmetry can 
be altered dramatically [j, HJ [l5|] • I n particular, at high 
q 2 , the sign of Afb is sensitive to the sign of the prod- 
uct of the Cg S and C*q Wilson coefficients. The value 
of Afb in B — ► K£ + £~ provides an important check 
on this measurement, as it is expected to result in zero 
asymmetry for all q 2 in the SM and many non-SM sce- 
narios. This condition can be violated in models in which 
new operators such as a neutral Higgs penguin contribute 
significantly |18j . However even in this case the resulting 
asymmetry is expected to be of order 0.01 or less in the 
B — * K£ + £~ mode for electron or muon final states |2Cl |. 
In addition to Afb, m B — ► K*£ + £~ the fraction of 
longitudinal polarization Fl of the K* can be measured 
from the angular distribution of its decay products. The 
value of Fl measured at low q 2 is sensitive to effects from 
new left-handed currents with complex phases different 
from the SM, resulting in Cf = -C 7 (SM), or effects 
from new right-handed currents in the photon penguin 
amplitude J21J. The predicted distributions of ApB{q 2 ) 
and Fi(g 2 ) are shown for the SM and for several non-SM 
scenarios in Figure [3 The non-SM scenarios correspond 
to those studied in Refs. 0, 0, H2 . 

Finally, the lepton flavor-violating decays B —* 
K^*>e (i^ can only occur at rates far below current ex- 
perimental sensitivities in the context of the SM with 



neutrino mixing. Observation of these decays would 
therefore be an indication of contributions beyond the 
SM. For example, such decays are allowed in leptoquark 
models Ki- 



ll. DETECTOR AND DATASET 

The results presented here are based on data collected 
with the BABAR detector at the PEP-II asymmetric e + e~ 
collider located at the Stanford Linear Accelerator Cen- 
ter. The datasct comprises 229 million BB pairs, cor- 
responding to an integrated luminosity of 208 fb" 1 col- 
lected on the Y(4S) resonance at a center-of-mass energy 
of y/s = 10.58 GeV. An additional 12.1 ftT 1 of data col- 
lected at energies 40 MeV below the nominal on-peak 
energy is used to study continuum backgrounds arising 
from pair production of u, d, s, and c quarks. 

The BABAR detector is described in detail in Ref. 123 • 
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FIG. 2: Simulated distribution of (a) Afb and (b) Fl for 
the decay B — > K*£ + £~. The points represent the distri- 
butions assuming the SM (solid lines), Cf = -C 7 (SM) 
(dotted lines), CfCIo = -CgCio(SM) (dashed lines), and 
Cf,CfC1o = -C 7 (SM),-C 9 Cio(SM) (dot-dashed lines) 
generated using the form factor model of 123 . In the case of 
Fl, the two solutions with Cg S Cio — -CqC w (SM) are not 
displayed; they are nearly identical to the two shown. 



The measurements described in this paper rely primarily 
on the charged particle tracking and identification prop- 
erties of the detector. Tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) in a 1.5-T magnetic field produced by a supercon- 
ducting magnet. Low momentum charged hadrons are 
identified by the ionization loss (dE/dx) measured in the 
SVT and DCH. and higher momentum hadrons by a ring- 
imaging detector of internally reflected Chercnkov light 
(DIRC). A CsI(Tl) electromagnetic calorimeter (EMC) 
provides identification of electrons, and detection of pho- 
tons. The steel in the instrumented flux return (IFR) 
of the superconducting coil is interleaved with resistive 
plate chambers, providing identification of muons and 
neutral hadrons. 



III. EVENT SELECTION 

We reconstruct signal candidates in eight final states: 
B+ -s. K+l+e-, B° -► K°£+£-, B° -> K* £ + £~, B+ -> 
K* + £ + £~, where K*° -> K + n~, K* + -> K°tt + , K° -> 
7r + 7r _ , and £ is either an e or fx. Throughout this paper, 
charge-conjugate modes are implied. 

Electrons are required to have momentum above 0.3 
GeV/c and are identified using a likelihood ratio combin- 
ing information from the EMC, DIRC, and DCH. Pho- 
tons that lie in a small angular region around the elec- 
tron direction and have E > 30 MeV are combined with 
electron candidates in order to recover bremsstrahlung 
energy. We suppress backgrounds due to photon con- 
versions in the B — > Ke + e~ channels by removing e + e~ 
pairs with invariant mass less than 0.03 GeV/c 2 . As there 
is a significant contribution to the B — > K*e + e~ channels 
from the pole at low dielectron mass, we preserve accep- 
tance by vetoing conversions in these channels only if the 
conversion radius is outside the inner radius of the beam 
pipe. Muons with momentum p > 0.7 GeV/c are identi- 
fied with a neural network algorithm using information 
from the IFR and the EMC. 

The performance of the lepton identification algo- 
rithms is evaluated using high-statistics data control 
samples. The electron efficiency is determined from sam- 
ples of e + e~ — > e + e~7 events to be approximately 91% 
over the momentum range considered in this analysis; 
the pion misidentification probability is < 0.15%, evalu- 
ated using control samples of pions from r and K® de- 
cays. The muon efficiency is approximately 70%, deter- 
mined from a sample of e + e~ — + fi + fi~"/ decays; the pion 
misidentification probability is of order 2 — 3%, as de- 
termined from r decays. These samples are used to cor- 
rect for any discrepancies between data and simulation as 
a function of momentum, polar angle, azimuthal angle, 
charge, and run period. 

Charged kaons are selected by requiring the Cherenkov 
angle measured in the DIRC and the track dE/dx to be 
consistent with the kaon hypothesis; charged pions are se- 
lected by requiring these measurements to be inconsistent 



with the kaon hypothesis. Kg candidates are constructed 
from two oppositely charged tracks having an invariant 
mass in the range 488.7 < m^n < 507.3 MeV/c 2 , a com- 
mon vertex displaced from the primary vertex by at least 
1 mm, and a vertex fit x 2 probability greater than 0.001. 
The K s mass range corresponds to a window of approxi- 
mately 3d about the nominal K s mass. Modes that con- 
tain a K* are required to have a charged K or Kg which, 
when combined with a charged pion, yields an invariant 
mass in the range 0.7 < m^^ < 1.1 GeV/c 2 . 

The performance of the charged hadron selection is 
evaluated using control samples of kaons and pions from 
the decay D° —> K~% + , where the D° is selected from 
the decay of a D*. The kaon efficiency is determined to 
be 80 — 97% over the kinematic range relevant to this 
analysis. The pion misidentification probability is < 3% 
for momenta less than 3GeV/c, and increases to ~ 10% 
at 5 GeV/c. As with the leptons, these samples are used 
to correct for any discrepancies between the hadron ID 
performance in data and simulation. 

Correctly reconstructed B decays will peak in two 
kinematic variables, tues and AE. For a candidate 
system of B daughter particles with total momentum 
p B in the laboratory frame and energy E* B in the 
T{AS) center-of-mass (CM) frame, we define tues — 
^{s/2 + p Q -p B ) 2 /E 2 -p 2 B and AE = E* B - ^i/2, 
where E and po are the energy and momentum of the 
T(4S) in the laboratory frame, and y/s is the total CM 
energy of the e + e~ beams. For signal events, the tties 
distribution peaks at the B meson mass with resolution 
cr w 2.5 MeV/c 2 . The AE distribution peaks near zero, 
with a typical width a ss 18 MeV in the muon channels, 
and a ps 22 MeV in the electron channels. 

B candidates are selected if the reconstructed tues 
and AE are in the ranges 5.00 < ities < 5.29 GeV/c 2 
and -0.50 < AE < 0.50 GcV. The signal is extracted 
by performing a multidimensional, unbinned maximum- 
likelihood fit in the region 5.20 < m E s < 5.29 GeV/c 2 
and -0.25 < AE < 0.25 GeV, which contains 100% 
of the signal candidates that pass all other selection re- 
quirements. This region remains blind to our inspection 
until all selection criteria are established. The events 
in the sideband with 5.00 < m E s < 5.20 GeV/c 2 , or 
-0.50 < AE < -0.25 GcV, or 0.25 < AE < 0.50 
GeV are used to study the properties of the combina- 
torial background. 

For the measurements of the partial branching frac- 
tions, Afb, and K* polarization, wc subdivide the sam- 
ple into two regions of dilepton invariant mass. The first 
is the region above the pole and below the J/tp reso- 
nance, 0.1 < q 2 < 8.41GcV 2 /c 4 ; the second is the re- 
gion q 2 > 10.24 GeV 2 /c 4 , above the J/tp resonance. The 
ip(2S) resonance is explicitly excluded from this upper 
region as described in further detail in Section liV Bl The 
lower bound of 0.1 GeV /c 4 in the first region is chosen to 
remove effects from the photon pole in the B — ► K*e + e~ 
channel. The forward-backward asymmetry is extracted 
in each of these q 2 regions from the distribution of cos 9* , 
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FIG. 3: Predicted distributions of (a) Afb{ci ) and (b) Fh{q ) in B — > K* £ + £ for the two regions of q considered. The lines 



off 



solutions with Co 



represent the predictions of the SM (solid lines), Cy = -C 7 (SM) (dotted lines), CgCfo = -C 9 Cio(SM) (dashed lines), and 
Cr(SM), —CgCio(SM) (dot-dashed lines) with the form factor model of Ref. I25I . In the case of Fl, the two 



ty^cff 

I-'IO — 



-CgCio(SM) are not displayed; they are nearly identical to the two shown. 



which wc define as the cosine of the angle between the 
£~ {l + ) and the B (B) meson, measured in the dilcp- 
ton rest frame. We do not measure Afb in the mode 
B° — > Kg£ + £~, in which the flavor of the B meson can- 
not be directly inferred from the K®. The K* polariza- 
tion is similarly derived from the distribution of cosOk, 
defined as the cosine of the angle between the K and the 
B meson, measured in the K* rest frame. The predicted 
distributions of Afb and Fl integrated over these two 
q 2 ranges are shown in Figure |2l for both the SM and 
non-SM scenarios. 



IV. BACKGROUND SOURCES 

A. Combinatorial backgrounds 

Combinatorial backgrounds arise either from the con- 
tinuum, in which a (u, d, s, or c) quark pair is produced, 
or from BB events in which the decay products of the 
two S's are mis- reconstructed as a signal candidate. We 
use the following variables computed in the CM frame to 
reject continuum backgrounds: (1) the ratio of second to 
zeroth Fox- Wolfram moments [26j , (2) the angle between 
the thrust axis of the B and the remaining particles in 
the event, Othrust, (3) the production angle Bb of the 
B candidate with respect to the beam axis, and (4) the 
invariant mass of the kaon-lepton pair with the charge 
combination expected from a semilcptonic D decay. The 
first three variables take advantage of the characteris- 
tic jet-like event shape of continuum backgrounds, ver- 
sus the more spherical event shape of BB events. The 
fourth variable is useful for rejecting cc events. These 
frequently occur through decays such as D — > K£v, re- 
sulting in a kaon-lepton invariant mass which peaks be- 



low that of the D; for signal events the kaon-lepton mass 
is broadly distributed up to approximately the B mass. 
These four variables are combined into a linear Fisher 
discriminant |27| . which is optimized using samples of 
simulated signal events and off-resonance data. A sep- 
arate Fisher discriminant is used for each of the decay 
modes considered in this analysis. 

Combinatorial BB backgrounds are dominated by 
events with two semilcptonic B — » Xiv decays. We dis- 
criminate against these events by constructing a likeli- 
hood ratio composed of (1) the vertex probability of the 
dilepton pair, (2) the vertex probability of the B candi- 
date, (3) the angle 0b as in the Fisher discriminant, and 
(4) the total missing energy in the event E m i SS . Events 
with two semilcptonic decays will contain at least two 
neutrinos; therefore the E m i SS variable is particularly ef- 
fective at rejecting these backgrounds. The probability 
distribution functions (PDFs) used in the likelihood arc 
derived by fitting simulated signal events and simulated 
BB events in which the signal decays are removed. We 
derive a separate likelihood parameterization for each de- 
cay mode. 

We select those events that pass an optimal Fisher and 
BB likelihood requirement, based on the figure of merit 
S/y/S + B for the expected number of signal events S 
and background events B. The selection is optimized si- 
multaneously for the Fisher and likelihood, and is derived 
separately for each decay mode. 



B. Peaking backgrounds 

Backgrounds that peak in the tues an d A-E vari- 
ables in the same manner as the signal are either ve- 
toed, or their rate is estimated from simulated data or 
control samples. The largest sources of peaking back- 
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FIG. 4: Charmonium veto regions (a) in the B + — > if + e + e _ channel. The points are simulated J ftp and ip{2S) events, with 
abundance equal to the mean number expected in 208 fb _1 . The projections onto (b) m e + e - and (c) AE are shown at right, 
indicating the high density of points at [m l + l - , AE) — (m^,0.0). The vertical band corresponds to events where the J/ip 
(ip(2S)) and K + come from different B decays. For AE < it also includes events with mis-reconstructed B — > J/ipK* , 
B — ► ip(2S)K* , and non-resonant charmonium decays. The slanted band corresponds to events with mis-measured lepton track 
momentum. 



grounds are B decays to charmonium: B — ► J/ipK^*' 
and B -> tp(2S)K^, where the J/ip or ^(2S) decays 
to a £ + ^~ pair. We therefore remove events in which 
the dilepton invariant mass is consistent with a J/ip or 
ip(2S), either with or without bremsstrahhmg recovery 
in the electron channels. In cases where the lepton mo- 
mentum is mis-measured, or the brcmsstrahlung recov- 
ery algorithm fails to find a radiated photon, the dilepton 
mass will be shifted from the charmonium mass. In addi- 
tion, the measured AE will be shifted away from zero in 
a correlated manner. We account for this by construct- 
ing a two-dimensional veto region in the mg+g- vs. AE 
plane as shown in Figure 0J the simulated points plotted 
demonstrate the expected background rejection. Within 
the veto region in data we find approximately 13700 J/ip 
events and 1000 ip(2S) events summed over all decay 
modes. These provide a high-statistics control sample 
useful for evaluating systematic uncertainties and selec- 
tion efficiencies. The residual charmonium background 
after applying the veto is estimated from simulation to 
be between 0.0 and 1.6 events per decay mode. 



Due to the 2-3% probability for misidentifying pions as 
muons, the B — * K^*' n + [i~ channels also receive a signif- 
icant peaking background contribution from hadronic B 
decays. The largest of these are B~ — ► D°tt~ where 



D° 
where D + 



D° 



K 



and B a 



D+tt- 



K*°7r + . These are suppressed by removing 
events in which the K^*'fi invariant mass lies in the range 
1.84 < m K {,) < 1.90 GeV/c 2 . The remaining hadronic 
backgrounds come from charmless decays such as B —> 



ifWvr+Tr-, B -> K^K+tt-, and B -► K^K+K'. We 
measure the peaking background from these processes us- 
ing data control samples of B — » K^hfi events. These 
samples are selected with the same requirements as sig- 
nal events, except hadron identification is required for 
the hadron candidate h in place of muon identification. 
This yields samples of predominantly hadronic B decays. 
We then weight each event by the muon misidentification 
rate for the hadron divided by its hadron identification 
efficiency. The hadronic peaking background is then ex- 
tracted by a fit to the ities distribution of these weighted 
events. This results in a total hadronic peaking back- 
ground measurement of 0.4 - 2.3 events per muon decay 
channel. These backgrounds are suppressed by a factor 
of approximately 400 in the B — > K(*^e + e~ channels due 
to the much lower probability of misidentifying pions as 
electrons. 

There is an additional contribution to the peaking 
backgrounds in the electron channels from rare two-body 
decays. These include B — > K*j with the 7 converting 
to an e + e~ pair in the detector, and B — > K^ir or 
B — > K(*'r), where the ir° or r\ undergoes a Dalitz decay 
to e + e~7. These backgrounds are estimated from sim- 
ulation to contribute 0.0 - 1.4 events per electron decay 
channel. 

The sum of peaking backgrounds from all sources is 
summarized in Table [I] As a function of q 2 , all of 
the backgrounds from K*^ and K^tt are localized in 
the region 0.0 < q 2 < 0.1 GeV /c 4 . Backgrounds from 
J/ip and K^rj populate the region 0.1 < q 2 < 8.41 



TABLE I: Mean expected peaking backgrounds in 208 fb -1 , 
for the individual K^*'£ + £~ decay modes after applying all 
selection requirements. 





Allg" 


0.1 < q 1 < 8.41 


q l > 10.24 


Mode 




(GeV 2 /c 4 ) 


(GeV 2 /c 4 ) 


K+e+e- 


0.7 ±0.2 


0.6 ±0.2 


0.1 ±0.1 


K + ii + pT 


2.3 ±0.5 


1.4 ±0.4 


0.9 ±0.1 


K%e+e- 


0.01 ±0.01 


0.01 ±0.01 


0.0 


K° sf i+n- 


0.4 ±0.1 


0.3 ±0.1 


0.1 ±0.04 


K*°e+e- 


3.0 ±0.6 


1.0 ±0.5 


0.6 ±0.2 


ff*W 


1.4 ±0.8 


0.5 ±0.3 


0.2 ±0.1 


K*+e+e- 


0.9 ±0.2 


0.2 ±0.2 


0.2 ±0.1 


K*+n+ii- 


0.6 ±0.3 


0.2 ±0.1 


0.2 ±0.1 



GeV /c 4 , while the ip{2S) backgrounds contribute only 
to the region q 2 > 10.24 GeV 2 /c 4 . The hadronic back- 
grounds occupy both the 0.1 < q 2 < 8.41 GeV 2 /c 4 and 
q 2 > 10.24 GeV 2 /c 4 regions. 



V. YIELD EXTRACTION PROCEDURE 

We extract the signal yield and angular distributions 
using a multidimensional unbinned maximum likelihood 
fit. For B — > K£ + £~ , the total branching fraction is ob- 
tained from a two-dimensional fit to tues an d AE. In 
the B — > K*£ + £r modes, we add the reconstructed K* 
mass as a third fit variable. The signal shapes are pa- 
rameterized in both tties and AE by a Gaussian function 
plus a radiative tail described by an exponential power 
function. This takes the form 



(x-x) 2 \ 



f(x)<xi eXP( , 2. 

' A x (B - ^ 



; (x — x)/<7 > a 

)"" ; {x-x)/a<a ' 



where A = (^) n x exp(-|a| 2 /2) and B = A - |a|. The 
variables x and a are the Gaussian peak and width, and 
a and n are the point at which the function transitions to 
the power function and the exponent of the power func- 
tion, respectively. The toes shape parameters x, a, a, 
and n are assumed to have a AE dependence of the form 
c o ± c-2{AE) 2 , determined empirically from simulation. 
The mean and width are fixed to the values derived by 
fitting the control sample of vetoed charmonium events. 
All other signal shape parameters are fixed to the val- 
ues obtained from fits to simulated signal events. In the 
B — > K*£ + £~ mode, the mass of the K* is parameterized 
with a relativistic Breit-Wigner line shape. 

The background is modeled as a sum of terms de- 
scribing (1) combinatorial background; (2) peaking back- 
ground; (3) cross-feed backgrounds; and, (4) in the 
B — > K*£ + £~ modes, backgrounds that peak in rriKir 



at the K* mass but not in rriEs and AE. The com- 
binatorial background is described by a product of an 
empirically derived threshold function in rriEs- a linear 
term in AE, and the product of \/itiktt — m-K — "m-ir and 
a quadratic function of toa'tt for the K* modes. The form 
of the threshold function used to describe the background 
in rriEs is f{x) oc xy/l — x 2 cxp [— £(1 — x 2 )], where £ is 
a fit parameter and x = toes/-Ev The peaking back- 
ground component has the same shape as the signal, with 
normalization fixed to the estimates of the mean peak- 
ing backgrounds (Table QJ. The cross-feed component 
has a floating normalization to describe (a) background 
in B -> K(+£- (B -► K*£+£-) from B -> K*l + t~ 
(B — ► K*tt£ + £~) events with a lost pion, and (b) back- 
ground in B — > K*£ + £~ from B — > K£ + £~ events with a 
randomly added pion. The backgrounds that peak only 
in mK-K & re described by the signal shape in juk n an d the 
combinatorial background shape in tties an d AE. The 
yield of this term is fixed to (5 ± 5)% of the total combi- 
natorial background, as determined from simulation. As 
the shape parameters for term (1) and the normalizations 
for terms (1) and (3) are all free parameters of the fit, 
much of the background uncertainty propagates into the 
statistical uncertainty in the signal yield obtained from 
the fit. 

The CP asymmetry is also extracted from the fit in 
the B+ -> K+£+£- and B -► K*£+£~ channels, where 
the flavor of the b quark can be inferred from the charge 
of the final state K^*> hadron. As this cannot be done 
in the case of B° — ► K®£ + £~, we do not measure the 
CP asymmetry in that mode. The possibility of a non- 
zero CP asymmetry in the combinatorial background is 
accounted for by allowing its value to float in the fit. 
The CP asymmetry of the peaking background is fixed 
to the value expected from the relative composition of 
background sources. 

The partial branching fractions are measured by re- 
peating the fit with the sample partitioned into q 2 bins. 
The signal efficiencies and peaking backgrounds are re- 
computed for each region of q 2 . To determine the 
forward-backward asymmetry and K* polarization in 
bins of q 2 , we also utilize fits to the cos 8* and cos 0k an- 
gular distributions. We follow the treatment of Ref. |21| 
to parameterize the angular distributions for signal. The 
signal shape in cos 9k is described by an underlying dif- 
ferential distribution which depends on the fraction of 
longitudinal polarization Fl as 



1 dT = 3 F c g2 
rdcos6 K 2 LC ° S 



IK 



-(l-F L )(l-cos 2 ^). 



The underlying differential rate for signal in cos 9* is 
then described in terms of F^ and the forward-backward 
asymmetry term Aeb which enters linearly in cos 9* : 



1 dT 



r d cos ( 



-F L {\- cos 2 9* 



±-(l-F L )(l±cos 2 r)± 
Aeb cos 8* . 
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In the B + — > K + £ + £ mode, the most general distribu- 
tion for cos 9* with non-zero Aps is given by: 



1 dT 



1 



T d cos ( 



-(1 -F s )(l - cos 2 9*) + -F s + A FB cos9 



where Fs is the relative contribution from scalar and 
pseudoscalar penguin amplitudes, and Afb arises from 
the interference of vector and scalar amplitudes J22| ■ In 
the Standard Model, both Fs and Afb are expected to 
be negligibly small; their measurement is therefore a null 
test sensitive to new physics from scalar or pseudoscalar 
penguin processes. 

The true angular distributions are altered by detec- 
tor acceptance and efficiency effects. We account for 
this by multiplying the underlying distributions with ef- 
ficiency functions e(cos9*) and c{cos9k) described by a 
non-parametric histogram PDF obtained from signal sim- 
ulations. 

The combinatorial background shapes in cos9* and 
cos 9 k are described by a histogram PDF drawn from 
control samples in the rriES and AE sidebands. The an- 
gular distribution of the peaking backgrounds are fixed 
in the fit. Additional components describing the angular 
distribution of cross-feed events and of mis-reconstructed 
signal events are included as histogram PDFs derived 
from simulated samples. 

In the B — > K*£ + £~ modes we first perform a four- 
dimensional fit to rriES, AE, rriKir, and cos9k to ob- 
tain Fl. Due to limited statistical sensitivity of Fl to 
the cos 9* distribution, Fl is fixed to the value mea- 
sured from the cos 9k distribution in order to measure 
Afb from a fit to niEs, AE. ttlkti, and cos 9*. In the 
B + — > K + £ + £~ modes, Afb and Fs are simultaneously 
extracted directly from a three-dimensional fit to rriESi 
AE, and cos6>*. 



VI. SYSTEMATIC UNCERTAINTIES 
A. Branching fractions 

In evaluating systematic uncertainties in the branch- 
ing fractions, we consider both errors that affect the sig- 
nal efficiency estimate, and errors arising from the max- 
imum likelihood fit. Sources of uncertainties that affect 
the efficiency are: charged-particle tracking (0.8% per 
lepton, 1.4% per charged hadron), charged-particle iden- 
tification (0.5% per electron pair, 1.3% per muon pair, 
0.2% per pion, 0.6% per kaon), the continuum back- 
ground suppression selection (0.3%-2.2% depending on 
the mode), the BB background suppression selection 
(0.6%-2.1%), K° selection (0.9%), and signal simulation 
statistics (0.4%-0.7%). The estimated number of BB 
events in our data sample has an uncertainty of 1.1%. 
We use the high-statistics sample of events that fail the 
charmonium veto to bound the systematic uncertainties 
associated with the continuum suppression Fisher dis- 
criminant, the BB likelihood suppression selection, and 



charged particle identification. The Fisher discriminant 
and BB likelihood ratio for B + — > K + e + e~ are illus- 
trated in Figure for data and simulation in the J/ip 
control sample. An additional systematic uncertainty in 
the efficiency results from the choice of form factor model, 
which alters the q 2 distribution of the signal. We take this 
uncertainty to be the maximum efficiency variation ob- 
tained from a set of recent models 0, Is Hoi l24l l25j ; the 
uncertainty is computed separately for each mode and 
varies in size from 1.1% to 8.3%. 

Systematic uncertainties on the signal yields obtained 
from the maximum-likelihood fit arise from three sources: 
uncertainties in the parameters describing the signal 
shapes, uncertainties in the combinatorial background 
shape, and uncertainties in the peaking backgrounds. 
The uncertainties in the means and widths of the signal 
shapes arc obtained by comparing data and simulated 
data in B — > J/ipK^*' control samples. For modes with 
electrons, we also vary the fraction of signal events in the 
tail of the AE distribution by varying the exponent n in 
the exponential power function. Signal shape uncertain- 
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FIG. 5: Distribution of (a) the Fisher discriminant and (b) 
the BB likelihood ratio for B + — > K + e + e~ events in the J/ip 
veto sample. The points are data; the gray bands are simu- 
lated events, with a simulation uncertainty given by the band 
height. The dark gray portion represents the uncertainty due 
to simulation statistics, while the additional uncertainty due 
to the B — » J/ipK + branching fraction is represented by the 
light gray band height. Events to the right of the vertical line 
are selected. 
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ties are typically 2-4% of the signal yield. To evaluate 
the uncertainty due to the background shape, we reeval- 
uate the fit yields with three different parameterizations: 
(1) an exponential shape for AE, (2) a quadratic shape 
for AE, and (3) an rriEs background shape parameter 
£ which is linearly correlated with AE. In modes with 
a K* , we also vary the yield of the background com- 
ponent which peaks in ttiktv but not in rriES or AE by 
100% of itself. The induced uncertainty in the signal yield 
due to the background shape is 4-6% for B — ► K£ + £~ 
modes and increases to 8-12% for B — > K*£ + £~ modes, 
where the backgrounds are generally larger. Uncertain- 
ties in the peaking background induce an uncertainty in 
the signal yields of 2-5%; this is obtained by varying 
the expected peaking background yield within its ±lcr 
uncertainties. The total systematic uncertainty in the 
fitted signal yield induces a systematic uncertainty ABm 
in the measured branching fraction; this uncertainty is 
shown for each of the branching fraction fits in Tables |n] 
and EH 



B. CP asymmetry 

The systematic uncertainties in the measurement of 
Acp include errors due both to detector efficiency ef- 
fects and to the asymmetry in the peaking background 
component. The error associated with the detector effi- 
ciency is obtained by comparing the value of Aqp mea- 
sured in the charmonium control samples with the ex- 
pected value of zero; agreement with zero is obtained 
with a precision of 1.2% for B+ -> K+£+£- and 2.1% 
for B — ► K*£ + £~ . The uncertainty due to the peak- 
ing background is evaluated by varying the expected CP 
asymmetry of the peaking backgrounds within their un- 
certainties. The possible CP asymmetry in the charmo- 
nium and B —> A"* 7 peaking backgrounds is highly con- 
strained from previous measurements; any asymmetry in 
the Dalitz decays is suppressed by their relatively small 
contribution to the peaking background. In contrast, the 
hadronic peaking background in the muon modes could 
exhibit a significant CP asymmetry; this is measured di- 
rectly from the asymmetry of the hadronic control sample 
described in Section flV Bl with an uncertainty dominated 
by the statistics of the sample. This induces an uncer- 
tainty in the measured Acp of 1% for B + — > K + £ + £~ 
and 2% for B — > K*£ + £~ . Other systematic uncertain- 
ties induced by the fitting procedure, as computed above 
for the branching fraction measurements, are found to be 
negligible. 



C. Angular distributions 

Systematic uncertainties related to the angular distri- 
butions of the efficiency arc estimated by comparing the 
values of Afb, Fg, and Fl measured in the relevant char- 
monium control samples with their expected values. For 



B — + J/ipK and B — > J/ipK* we measure an App con- 
sistent with zero and with a precision of 0.01 and 0.02, 
respectively. For B — ► J/tpK*, we measure Fl to be 
consistent with the previous BABAR measurement |28j . 
with a precision of 0.05. For B — ► J/ipK we measure F$ 
consistent with zero and with a precision of 0.03. 

Further systematic uncertainties arc evaluated by re- 
peating the fit with alternative shapes assumed for 
the background components: (1) the shape of mis- 
reconstructed signal events is fixed instead to the shape 
of correctly reconstructed signal, (2) the combinatorial 
background shape is drawn from alternative ranges of 
rriES an d AE, and from the sample of events that fail 
the BB likelihood selection, and (3) the angular distri- 
butions of the peaking backgrounds arc varied within 
their statistical uncertainties. Systematic uncertainties 
from backgrounds induce uncertainties in Fl and App 
of 0.05-0.18, depending on the relative amount of back- 
ground, and are the largest systematic uncertainty. F$ is 
more sensitive to the background shape, with an induced 
systematic uncertainty of 0.45. 

In the fit to cos 6* in the B — > K*£ + £~ decay modes, 
the value of Fl is fixed to the result obtained from the fit 
to the cos 9k distribution. This introduces an additional 
parametric uncertainty of 0.01 on the measured value 
of Afb , which we evaluate by varying Fl within the 
uncertainty of the measurement. 



VII. RESULTS 

A. Branching fractions 

We first perform the fit separately for each of the 
eight decay modes to extract the branching fractions inte- 
grated over all q 2 . In the branching fraction fits, the effi- 
ciency is defined such that the total branching fraction in- 
cludes the estimated signal that is lost due to the charmo- 
nium vetos. The results for the individual decay modes 
are shown in Table [H] We then perform a combined fit 
to the appropriate combinations of modes to extract the 
B — ► K£ + £~ and B — > K*£ + £~ branching fractions. We 
combine charged and neutral modes by constraining the 
total width ratio T(B°)/T(B + ) to the world average ra- 
tio of lifetimes t{B+)/t{B°) = 1.071 ±0.009 Hj. In the 
B — > K*£ + £~ mode, we add the additional constraint 
T(B -s. K*i_i+fi-)/T(B -► K*e+e-) = 0.75 to account 
for the enhancement due to the pole at low q 2 in the 
electron channel [j]. The final branching fractions are 
expressed in terms of the B° — > A"'*' (i + \i~ channels. 
With these constraints, we find the lepton-fLavor aver- 
aged, i3-charge averaged branching fractions 



B(B -> K£+£-) = (0.34 ± 0.07 ± 0.02) x 10" 



B{B -> K*t 



(0.78±g;^±0.11) x 10~ 6 , 



where the first error is statistical and the second system- 
atic. The projections of the data overlaycd with the com- 
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bined fit results are shown in Fi gures [S] an d [7| The signal 
significance is computed as ^/2Aln(£), where Aln(£) 
is the difference between the likelihood of the best fit 
and that of the null signal hypothesis. Systematic un- 
certainties are incorporated in the significance estimate 
by simultaneously applying all variations that result in a 
lower signal yield before computing the change in likeli- 
hood. The significance of the signal including statistical 
and systematic uncertainties is 6.6 standard deviations 
for the B — > Kl + t~ mode and 5.7 standard deviations 
for the B — ► K*£ + £~ mode. The secondary peak in the 
AE sideband of B — > K£ + £~ results from the fit com- 
ponent describing events with a lost pion, either from 
B — > K*£ + £~ or from events in which a b — ► s£ + £r de- 
cay results in a Kn£ + £~ final state without proceeding 
through an intermediate K* resonance. The normaliza- 
tion and mean AE of this component are free parame- 
ters in the fit. Examination of these events shows that 
the addition of a charged or neutral pion results in a 
B -> K*£+£r or B -> Ktt£+£- signal candidate. Using 
simulated signal decays, we find the effect of these events 
on the B — ► K£ + £~ signal yield is negligible. 

We further perform a set of combined fits with the 
sample partitioned into final states containing muons and 
electrons, and into charged and neutral final states, mod- 
ifying the constraints as appropriate. The results from 
all such fits are summarized in Table ITTT1 



TABLE II: Results from fits to the individual K^t+f de- 
cay modes for all q . The columns from left are: decay mode, 
fitted signal yield, signal efficiency, relative uncertainty on 
the branching fraction due to the systematic error on the effi- 
ciency estimate, systematic error from the fit, and the result- 
ing branching fraction (with statistical and systematic errors). 
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If the pole region is removed by requiring q > 
0.1GeV 2 /c 4 , the constrained ratio between B —* 
K*/j, + fi~ and B — > K*e + e~ in the combined fit is modi- 
fied from 0.75 to 1. Repeating the combined fit with this 
modification, we obtain 

B(B -f K*£ + r) {q2>0 , 1GcV 2 /ci) = (0.73™±0.11)xHT e 

The results of the combined fits in the various subsam- 
ples with the pole region removed are shown in Table ITTT1 
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FIG. 6: Distributions of the fit variables in K£ + £~ data 
(points), compared with projections of the combined fit 
(curves): (a) toes distribution after requiring —0.11 < AE < 
0.05 GeV and (b) AE distribution after requiring |m,ES — 
tub] < 6.6 MeV/c 2 . The solid curve is the sum of all fit com- 
ponents, including signal; the dashed curve is the sum of all 
background components. 
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(points), compared with projections of the combined fit 
(curves): (a) tues after requiring —0.11 < AE < 0.05 GeV 
and 0.817 < m K n < 0.967 GeV/c 2 , (b) AE after requiring 
|m ES - m B \ < 6.6 MeV/c 2 , 0.817 < m K * < 0.967 GeV/c 2 , 
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-0.11 < AE < 0.05 GeV. The solid curve is the sum of all 
fit components, including signal; the dashed curve is the sum 
of all background components. 



We observe good agreement in the branching fraction ob- 
tained in all of the subsamples, both with and without 
the pole region included. The measured total rates are 
consistent with the range of Standard Model rates pre- 
dicted in Ref. 0. The B -> K£+£~ rate is significantly 
lower than the range given by Ref. |12j . 

From the separate fits to the muon and electron chan- 
nels integrated over all q 2 , we obtain the ratios 

R K = 1.06 ±0.48 ±0.08, 
R K * =0.91 ±0.45 ±0.06, 
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TABLE III: Results from fits to combined K^l+T decay 
modes for all q 2 . The columns from left are: decay mode 
combination, fitted signal yield, relative uncertainty on the 
branching fraction due to the systematic error on the effi- 
ciency estimate, systematic error on the branching fraction 
introduced by the systematic error on the fitted signal yield, 
and the resulting branching fraction (with statistical and sys- 
tematic errors). The constraints for each combined fit are 
described in the text. 



Mode 


Yield 
(events) 


AB cff 
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0.73+.g;g ± 0.21 
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0.78±gif ±0.11 


Pole excluded 
K*e + e~ 
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±5.2 


±0.11 


0.63±o;22±0.11 


K'fi+fi- 


20.7±fo 


±5.9 


±0.11 


0.881^0 ±0.12 


K*°£+£- 


34.8i 9 ° 3 4 


±5.0 


±0.10 


0-75±oJo±0.10 


K*+£+£~ 


9.511;? 


±6.6 


±0.19 


0.73tog±0.19 


K*£+£- 


44.3±iii 


±5.3 


±0.11 


0.73+£?8 ± 0.11 



consistent with the SM predictions of 1.00 and 0.75, re- 
spectively. If instead the pole region is excluded from the 
B — > K*£ + £~ channels, we find 



C. Partial branching fractions 

The partial branching fractions obtained from the fits 
to nigs, AE, and mx^ in two bins of q 2 are shown in 
Table llVl The results are generally consistent with the 
q 2 dependence predicted in recent Standard Model based 
form factor calculations (Figure |SJ). 



D. K* polarization 

The fit projections for the cos 9k distribution in bins 
of q 2 are shown in Figure IT21 of Appendix 1X1 The result- 
ing values for the fraction of longitudinal polarization 
Fl are listed in Table IIVI Combining all events with 



q 2 > O.lGeVVc 4 , we find 



F L (B -+ K*t 



'( 9 2 >0.1GeV 2 /c 4 ) 



0.63 



+0.18 
-0.19 



±0.05, 



where the first error is statistical, and the second system- 
atic. 

The measured values of Fl are consistent with the SM 
expectation in both q 2 ranges (Figure 0) and integrated 
over all q 2 > 0.1GeV 2 /c 4 . However, the large statistical 
uncertainties do not allow the determination of the sign 
of Cj from this measurement at present. 



E. Lepton forward-backward asymmetry 

The fit projections for the cos9* distribution in the 
B + — > K + £ + £~ mode are shown in Figure El of Ap- 
pendix^] Combining all events with q 2 > 0.1GeV 2 /c 4 , 
we find for the B + — ► K + £ + £~~ mode 



A FB (B + -» K+£ + £-) 



(<? 2 >o.iGov 2 /c 4 ) ~~ 0.15_ ; 23 ±0.08, 



^., ((? 2 >0.1 GeV 2 /c 4 ) = 1-40 ± 0.78 ± 0.10, 
where this ratio is expected to be 1 in the SM. 

B. CP asymmetry 

From the fit to the combined modes integrated over all 
q 2 , we find the direct CP asymmetries 

A CP (B+ -c !<+£+£-) = -0.07 ± 0.22 ± 0.02, 



A CP {B -> K*£+£-) = ±0.03 ± 0.23 ± 0.03, 

where the first error is statistical and the second system- 
atic. The measured values in both channels are consis- 
tent with the SM expectation of a negligible direct CP 
asymmetry. 



F S (B+ ^ K+£ + £-) {q2>01GcV2/ci) 



OMtoil ± 0.46, 



where the first errors are statistical, and the second sys- 
tematic. The correlation coefficient between these two 
meaurements is ±0.23. Both Afb and F$ are consistent 
with the SM prediction of zero. As a cross-check, we have 
also performed similar fits in the low and high q 2 regions 
for Apb , where due to limited statistics F$ must be fixed 
to zero; the resulting asymmetries are — 0.49+Qgg ± 0.18 
and 0.26+q' 2 ,! ± 0.03, respectively, which again are both 
consistent with zero asymmetry. 

The fit projections for the cos 9* distribution in the 
B — > K*£ + t~ mode arc shown in Figure El of Ap- 
pendix ^ and the resulting values of Afb listed in Ta- 
bic IIVI We find a large positive asymmetry in the high 
q 2 region, consistent with the SM expectation. This dis- 
favors new physics scenarios in which the product of the 
Cg S and C°o Wilson coefficients have the same magni- 
tude but opposite relative sign as in the SM, which would 
result in a large negative asymmetry at high q 2 (Fig- 
ure [SJ. 
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FIG. 8: Partial branching fractions in bins of q for (a) B — > K£ + £~ and (b) B — > K*£ + l~ , normalized to the total measured 
branching fraction. The points with error bars are data, the lines represent the central values of Standard Model predictions 
based on the form factor models of Refs. [24, 25] (solid lines), 11(1 (dashed lines), and 0,13 (dot-dashed lines). 



TABLE IV: Results from fits to the combined K^*'£ + £~ decay modes in bins of q 2 . The columns from left to right are: fitted 
q 2 range, partial branching fraction, longitudinal K* polarization Fl, and the lepton forward-backward asymmetry Afb- The 
first and second uncertainties are statistical and systematic, respectively. In B — > K£ + £~ , Afb is measured in the charged B 
decay modes only. The constraints for each combined fit are described in the text. The partial branching fractions are defined 
such that they include the estimated rate within the vetoed J/ip and ip(2S) resonance regions where appropriate. 
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0.1 -8.41 
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0.10l££t±0.01 
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-0A9t°o.ll ± 0.18 



>0.1 



0.34 



+0.07 



±0.02 



0.81+o;ei ± 0.46 



0.15^; 2 3 ±0.08 



For the low q 2 region and the region integrated over all 
q 2 > 0.1GeV 2 /c 4 , the Afb value corresponding to the 
maximum likelihood is positive, but is near the bound- 
ary at which a larger Aps will result in a negative, un- 
defined value for the extended likelihood function. For 
these maximally asymmetric cases the Ap b result is com- 
puted as a one-sided lower limit using a toy Monte Carlo 
method. For fixed values of Apb, we randomly gener- 
ate from the experimentally measured PDFs an ensem- 
ble of toy experiments, and find the value of Afb for 
which 5% of experiments in the ensemble have a maxim- 
ium likelihood fit resulting in a maximally positive Aps ■ 
The uncertainties in the other PDF parameters are ac- 
counted for by varying them randomly for each generated 



experiment in the ensemble according to normal distri- 
butions determined by the parameters' measured central 
values and uncertainties. We account for systematic un- 
certainties that do not correspond to continuous PDF pa- 
rameters, such as the choice of combinatorial background 
PDFs for cos6>*, by generating ensembles for each PDF 
variation and choosing that which results in the lowest 
lower limit. With this method, we find Aps > 0.19 at 
95% CL for the low q 2 region. Combining all events with 
q 2 > 0.1 GeV 2 /c 4 , we find for the B -> K*t + f mode at 
95% CL 

A FB (B - K*£ + n^ >0 .iGcvyc-) > 0-55. 
The corresponding fit projections shown in Figure [T4l are 
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FIG. 9: (a) Fi,(q ) and (b) AFB(q ) in B — > K*£ + £ . The points with error bars are data, with the arrow at low q in Afb 
indicating the 95% CL allowed region. The lines represent the predictions of the SM (solid lines), C7 = —Ci^SM) (dotted 
lines), CfCto = -CgCi {SM) (dashed lines), and Cf,CfCf$ = -C 7 (SM), -C 9 C W (SM) (dot-dashed lines) with the form 
factor model of Ref. l2al . In the case of Fl, the two solutions with Cg Cfo = — CgCio(SM) are not displayed; they are nearly 
identical to the two shown. 



produced by fixing the Afb of the signal component to 
its maximum physical value. 



F. Search for lepton flavor-violation 

We extract the signal yield in the B — ► Kefj, and 
B — * K*e/j, final states in a similar manner as the 
K(*)£+£- decays, with the particle identification require- 
ments modified to select e /j+ pairs. The signal effi- 
ciencies for these modes are determined from simulations 
where the B decays according to a simple three-particle 
phase space model. The results are shown in TablclVl As 
any physics that allows these decays will not necessarily 
affect the e+ \i~ and e _ /i+ states equally, we quote results 
for each charge state in addition to combined charge- 
averaged results. The projections of the data overlayed 
with the results of the combined fits are shown in Fig- 
ures ED and ^2 We find no evidence for a signal in any 
of these channels, and therefore set upper limits on these 
processes. For the combined lepton-charge averaged, B- 
charge averaged modes we find 

B(B-> Ken) <3.8x 1CT 8 , 



TABLE V: Results from fits to lepton flavor-violating decay 
modes. The columns from left are: decay mode, fitted signal 
yield, selection efficiency, relative uncertainty on the branch- 
ing fraction due to the systematic error on the efficiency esti- 
mate, systematic error on the branching fraction introduced 
by the systematic error on the fitted signal yield, and the 
90% C.L. limit on the branching fraction. The constraints for 
combined fits are described in the text. 

Mode Yield T(%) 23(10"*) g UL (10~ a )" 

9.1 
13 

9.1 
27 
53 
34 
58 

130 
99 

140 



K + e + n~ 


-3.5±?;i 


12.6 


-12.1+.^ ±2.3 


K+e-(i+ 


-0.8+iJ 


12.6 


-2.9+1* ±1.9 


K+efi 


-3.2+.?;? 


12.6 


— n.i±|:| =h 3.2 


K° ef j. 


oq+1.9 
^• y -1.3 


12.5 


-30± 13 ± 15 


K*°e+f,- 


1 1+3.6 
1 - 1 -2.1 


10.4 


7i f 3 ± 5 


K*°e-^+ 


1 1+3.5 
i - 1 -2.2 


10.4 


-7+ 14 ± 7 


K*°efi 


0.9+1° 


10.4 


el ;s± 9 


K*+e + n- 


0.4+ij 


10.0 


9+ 11 ± 22 


K* + e-fj,+ 


1 7+3.3 

i -'-2.0 


10.0 


-32± f s ± 15 


K*+efj, 


-0.2+1? 


10.0 


-it 59 ± 32 



Ke/j, 
K*en 


-4.9+. r j 

1.01" 


- -12.1±^±3.0 
48+. i 7 ± 11 


3.8 
51 



B(B -> K*e^) < 51 x 10" 8 , 

at 90% CL. These limits are significantly more stringent 
than those of previous searches |30t l31| . 



VIII. CONCLUSIONS 

We have measured the branching fractions, partial 
branching fractions, direct CP asymmetries, ratio of 



muons to electrons, fraction of longitudinal K* polariza- 
tion, and lepton forward-backward asymmetries in the 
rare FCNC decays B -> Kl + i~ and B -► K*i+i~. 

The branching fraction, Acp, Rk, and Fl results are 
all consistent with the Standard Model predictions for 
these decays. The values of Afb and the scalar con- 
tribution Fg measured in the B + — > K + £ + £~ channel 
are consistent with the expected value of zero. In the 
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1.2 5.22 5.24 5.26 5.28 
m ES (GeV/c 2 ) 



0.1 0.1 

AE(GeV) 



FIG. 10: Distributions of the fit variables in Ke^i data 
(points), compared with projections of the combined fit 
(curves): (a) toes distribution after requiring —0.11 < AE < 
0.05 GeV and (b) AE distribution after requiring |m,ES — 
tub] < 6.6 MeV/c 2 . The solid curve is the sum of all fit com- 
ponents, including signal; the dashed curve is the sum of all 
background components. 
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FIG. 11: Distributions of the fit variables in K*e[i data 
(points), compared with projections of the combined fit 
(curves): (a) tues after requiring —0.11 < A_E < 0.05 GeV 
and 0.817 < m Kn < 0.967 GeV/e 2 , (b) AE after requiring 
|m ES - m B \ < 6.6 MeV/c 2 , 0.817 < m K n < 0.967 GeV/c 2 , 
and (c) mKir after requiring |m,ES — ttib\ < 6.6 MeV/c 2 and 
-0.11 < AE < 0.05 GeV. The solid curve is the sum of all 
fit components, including signal; the dashed curve is the sum 
of all background components. 



B — > K*£ + £~ channel the large positive value of Afb 
at high q 2 is consistent with the SM and disfavors new 
physics scenarios in which the relative sign of the product 
of the Cg and Cio Wilson coefficients is opposite that of 
the SM. At low q 2 a positive value of Afb is also favored, 
with a 95% CL lower limit that is slightly above the SM 
prediction, as derived using the form factor models of 
Rcfs. USUI. 

In addition, we have obtained upper limits on the lep- 
ton flavor-violating decays B — ► Kefj, and B — > K*e/j, 
that are approximately one order of magnitude lower 
than those of previous searches. 

We note that the Belle collaboration has recently re- 
ported |32( a measurement of the integrated forward- 
backward asymmetries, finding Afb{B + — > K + £ + £~) = 
0.10 ±0.14 ±0.01 and A FB (B -> K*(+£') = 0.50 ± 
0.15 ±0.02. From a fit to the cosO* and q 2 distributions, 
they conclude that scenarios in which the product of Cg 
and Cio has the opposite sign as expected in the SM are 
disfavored, consistent with the results reported here. 

All of the measurements reported here are limited by 
statistical uncertainties, and can be improved with the 
addition of more data. 
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APPENDIX A: FITS TO ANGULAR 
DISTRIBUTIONS 

In this appendix we present plots of the cos 9k and 
cos 9* distributions in data, together with the projections 
of the combined fits used to extract Fl and Afb ■ Fig- 
ure El shows the fitted cos 9k distributions for each of 
the q 2 bins considered in this analysis. Figures IT^landlTll 
display the fitted cos 6* distributions for each of the q 2 



ranges for the B+ -> K+£+£~ and B -> K*£+£~ decay 
modes, respectively. For the fits to the cos 9* distribu- 
tions in the B —* K*£ + £~ mode, the K* polarization Fl 
is fixed to its measured value, as described in the text. 
The deviations from a smooth parabolic shape in the sig- 
nal component are the result of the efficiency and accep- 
tance corrections, which arc described by non-parametric 
histogram PDFs. 
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FIG. 12: Distributions of the fit variable cos 0k in B — > K*£ + £ data (points), compared with projections of the combined fit 
(curves) after requiring -0.11 < AE < 0.05 GeV, |m E s - m B \ < 6.6 MeV/c 2 , and 0.817 < m Kn < 0.967 GeV/c 2 . The solid 
curve is the sum of all fit components, the dashed curve is the sum of all background components, and the dot-dashed curve 



is the signal component, 
shown. 



The q 2 regions (a) 0.1 < q 2 < 8.41GeV 2 /c 4 , (b) q 2 > 10.24 GeV 2 /c 4 , and (c) q 2 > 0.1GeV 2 /c 4 are 
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FIG. 13: Distributions of the fit variable cos#* in B + —> K + l + l~ data (points), compared with projections of the combined 
fit (curves) after requiring —0.11 < AE < 0.05 GeV and |tties — itib\ < 6.6 MeV/e 2 . The solid curve is the sum of all fit 
components, the dashed curve is the sum of all background components, and the dot-dashed curve is the signal component. 
The q 2 regions (a) 0.1 < q 2 < 8.41 GeV 2 /c 4 , (b) q 2 > 10.24 GeV 2 /c 4 , and (c) q 2 > 0.1GeV 2 /c 4 are shown. The combined fits 
shown for (a) and (b) are performed by fixing Fs to zero. 




FIG. 14: Distributions of the fit variable cos (9* in _B — > K*£ + £ data (points), compared with projections of the combined fit 
(curves) after requiring -0.11 < AE < 0.05 GeV, \m ES - m B \ < 6.6 MeV/c 2 , and 0.817 < m Kn < 0.967 GeV/c 2 . The solid 
curve is the sum of all fit components, the dashed curve is the sum of all background components, and the dot-dashed curve 
is the signal component. The q 2 regions (a) 0.1 < q 2 < 8.41GeV 2 /c 4 , (b) q 2 > 10.24 GeV 2 /c 4 , and (c) q 2 > 0.1GeV 2 /c 4 are 
shown. The combined fits shown for (a) and (c) are performed by fixing Afb to its maximal physical value. 



